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Abstract 

We explicitly construct a diffeomoiphic pair (p(x),p~ l (x)) in terms of an appropriate 
quadric spline interpolating the prime series. These continuously differentiable functions are 
the smooth analogs of the prime series and the prime counting function, respectively, and 
contain the basic information about the specific behavior of the primes. We employ p~ x {x) 
to find approximate solutions of Diophantine equations over the primes and discuss how this 
function could eventually be used to analyze the von Koch estimate for the error in the prime 
number theorem which is known to be equivalent to the Riemann hypothesis. 



1 Introduction 

We shall use the following notation: N is the natural numbers set, p n or p{n) is the n-th prime, 
Pn — Pn+i — jt?„ — 1 is the number of composites in the interval (p n ,p n+ \), P is the prime numbers 
set, tz(x) is the prime counting function, Li(x) is the logarithmic integral , i.e., 



X 




The main objective of this paper is to show that the following pair of one-to-one mappings 

p(n):N->P, P~\q) :P^N 
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could be extended to a pair of diffeomorphisms over the real semi-axis (0, °°). 

Definition 1 The pair of functions f(x) G C^^(0,oo), g(x) G C^^(l,oo) is called prime number 
diffeomorphic if the following conditions hold 

f(n) = p n , «GN (2) 

f( n +l) = l(Pn+Pn+l), neN (3) 



2) 2 

= x VxG(l,oo), (4) 

g(f(x)) = x V^G(2,oo), (5) 

7t(x) = [g(x)\. (6) 

The diffeomorphisms f(x) and g(x) are called the prime curve and the prime counting curve, 

respectively. 

The function Kr(x) of Riemann-Von Mangoldt ([1], p. 34, (2), (3)) which can be expressed in 
terms of the zeros of the Riemann zeta function is the closest, among all known function, to 
the prime counting curve. However, the function Kr(x) is not invertible and cannot be used for 
the correspondence to the appropriate prime curve. It turns out that one can take the opposite 
way: construct an invertible interpolation of the prime series and then obtain from it a smooth 
counting curve. In this paper we describe such an invertible interpolant which is found among 
differentiable polynomial splines of minimal degree. 

These diffeomorphisms allows us to use in a natural way Fourier analysis as well as iterative 
methods for the solution of nonlinear problems in the case when it is necessary to account for the 
specific character of the non asymptotic behavior of the primes. As an example of the application 
of the above diffeomorphisms we consider in Sect. 3 an approximate method for the solution of 
Diophantine equations over P. 

We shall prove in Sect. 4 that the differentiable function p~ l (x) has the same asymptotic 
behavior as n(x) for x — > °°. This could be used for the analysis of the validity of the von Koch 
estimate in the form \p~ l (x) — hi(x)\/ ^\n(x) ~ const., which is known to be equivalent to the 
Riemann hypothesis (RH) [1,2]. 

2 Prime number diffeomorphisms based on a quadric spline 

Let us define the following functions 

a n( x ) = -2p n -i(x-n) 2 + (x-n) + p n , x>0, n = 2,3,..., 
a n( x ) = 2 Pn[x-n--) +(2p n +l) (x-n- -) + Pr ' r "" rl 



Their derivatives are 



da n (x) 



2 x 1 ' \ 2 2 



4p n _ l (n-x) + l, ^M = 4p n (x-n) + l. 



dx ! " dx 
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For any n = 2, 3, ... the functions are sewed together 



a n (n) 
da~ (x) 



— n + 



"n+1 

da~ +1 (x 



dx 
n + 



a n \ n ) = Pn, 

da+(x) 



x=n 
1 



dx 



dx 



at I n + 



da+(x) 



= 1. 



x=n 
1 



dx 



x—n+ 



iPn+Pn+l), 



= 2p n +l. 



(7) 
(8) 

(9) 
(10) 



X—II+ if 

Equations (7), (8), (9) and (10) define the following continuously differentiable quadric spline 



p(x) 



x+1, 0<x<j, 



a n (x), n — j<x<n, n = 2,3, 
a+(x), n<x<n + j, n = 2,3, 



(ID 



with first derivative 

dp(x) 
dx 



1, 0<x<§, 

4p n -l(n — x ) + 1) n — j<x<n, n = 2,3,.... 

4p n (x — + n<x<n + j, n = 2,3,... 



(12) 



Inverting the function a n (x) in the interval n — ^ < x < n and a+ (x) in the interval n < x < n + ^ 
gives the following inverse functions and their derivatives 

\-{%p n - i {p n -x) + \y 



b-{x) 

db~{x) 
dx 

dx 



n+- 



= n + 



4p n -i 

{Spn(x- Pn ) + 1)^-1 
Ap n 



(8p„_i(p„-x) + l) \ 

{sp n (x-p n ) + iy^. 



The functions b n (x), b^(x) and their derivatives are sewed together in a similar way like Eqs. (7), 
(8), (9) and (10): 



b n (Pn) 

db~{x) 



bn(Pn)=n, 

db+(x) 



dx „ 

Pn+Pn+l 



dx 



= 1, 



X=p n 



, + ( Pn + Pn+l\ 1 



db n+l (x) 



dx 



db+{x) 



_Pn+Pn+]_ 



dx 



x= Pn+Pn+\ 2p n + l 



(13) 
(14) 

(15) 

(16) 
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Finally Eqs. (13), (14), (15) and (16) define the continuously differentiable inverse spline 



x— 1, 



P-\X)- 

with first derivative 

dp- l (x) 
dx 



n + 



n + 



l-(8p„-i(p„-x)+lp 
(Sp„(x-p„)+l)?-l 



(&p„-l(pn-x) + l) ? 

(Sp n (x- Pn ) + 1)-K 



1 <X< §, 

^^<^<P„, n = 2,3,..., 
p n <-*<£«±!»±! 5 n = 2,3,..., 



1 <x< §, 

Ej^+EiL< x < Pnj n = 2,3,..., 
Pn < x <Pn±£n±L^ n = 2,3,... . 



Lemma 1 The derivatives of p{x) and p (x) satisfy the following inequalities 

l<fM£)<oo, *>0, 



— dx 



0< 



dp-\x) 
dx 



< 1, X>1. 



Proof: The above inequalities follow directly from the definitions (12) and (18). 

□ 

In the rest of this section we shall prove the following 
Theorem 1 

(i) The pair (p(x),p~ l (x)) is prime number diffeomorphic. 

(ii) The specific behavior of the prime and counting curves are traced by the invariants: 



1 = 



dp(x) 



dx 



dp 1 (x) 



■ 1 = sign 



- 1 = sign 



x=n 
2 



d p(x) 



dx 2 



d 2 p~\x) 



dx 



x=n—0 



n = 2,3,..., 



*=Pn 



sign 



d 2 p(x) 



dx 2 



sign 



dx 2 
d 2 p-\x) 



x=n+0 
2„-\ 



x=p n -0. 



dx 2 



X=p n +0 , 



Proof 
(i) 



(17) 



(18) 



(19) 
(20) 



n = 3,4,..., 

n = 3,4,... . 



(21) 
(22) 

(23) 



According to the definitions (11), (12), (17) and (18) we have the inclusions p(x) E 
CW(0,°o) and p~ l (x) e C^\l,o°). The validity of the interpolation conditions (2) and 
(3) follows from Eqs. (7) and (9). The mutual invertibility of p(x) and p~ l (x) follows 
from the fact that these functions are continuous and monotonically increasing (see (19) 
and (20)) and the conditions (4) and (5) can be checked directly. Equations (13) show that 
the functions 7t(x) and p~ l (x) are related by the identity (6). 
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(ii) Equations (21) follow from Eqs. (8) and (14). The discontinuity of the second derivatives 



~~ — 4p M _i, n — ^ < x < «, n = 2,3,... 

4/? M , n<x<n + i, n = 2,3,... 



d* 2 



and 



j2 — 1/ \ I 0, 0<x<2, 
" 7 ' U ' - / 4pn-i(%p n -i(p n -x) + l)-H\ ^Pjl< x < Pu , n = 2 ,3,..., 
-4^(8p„(x-p„) + l)- 3 / 2 , Pn <x<^±i ,n = 2,3,..., 



d* 2 



implies Eqs. (22) and (23). This completes the proof. 

□ 

Remark 1 The fact that the coefficients ofp(x) are integers as well as its invertibility follow from 
Eq. (3) in Definition 1 (i.e., from Eqs. (9)). 
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The comparative plot of the functions p~ l (x) and Kr(x) is shown on Fig. 1. The derivative 
dp~ x /dx is shown on Fig. 2. Its oscillating nature as well as the fact that it takes values between 
and 1 is obvious from this figure. 



1-0-. 




-i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — | — i — r 
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Figure 2: The derivative of p 1 (x) 



The diffeomorphisms p(x), p l (x) and their derivatives are realized in a Fortran90 program 
package called pp_ . f 90 which can be found at this URL [3]. 



3 Approximate solution of Diophantine equations 

Definition 2 Given a set of strictly increasing functions hi(x) G (0, °°), i = 1 , . . . n the system 

f 1 {x l ,...,x n ) = 0, (24) 

n 

f 2 {h{xi),. . . , h n (x n )) := £ sin 2 (7T^)) = 0, (25) 

(=l 
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where Eq. (24) is a Diophantine one is called real-Diophantine on the real semi axis (0,°°). 



The real-Diophantine systems allow us to find solutions of Diophantine equations in terms of real 
approximations of integer numbers. For this purpose one should apply numerical methods which 
work even when the derivative is degenerate at the solution (see, e.g., [4] and [5]). 
Let us write the system (24), (25) in a vector form as follows 

Fx = 0, (26) 

where 

Fx = f' T (x)fx, fx=[f 1 (x),f 2 (x)} T , /:fl/Cl^l m , xGK", n > m, 

f'(x) is the Jacobi matrix and Df is an open convex domain in R w . 

Here we shall quote an autoregularized version of the Gauss-Newton method [6, 7] as one of 
the possible methods for the solution of Eq. (26): 

x° e D f , e > 0, [f' T (x k )f'(x k ) + £*/) (x k+l -x k ) = -Fx k , k = 1,2, . . . . (27) 

£ k= \ (^/?k+ 4c Pk-?k^ , 

T*=||/'V)/V)|k P*=II^IU C=^±^, 

Po 

where 1 1 • | ]«, is the uniform vector or matrix norm, and the SVD method [8] is assumed for the 
solution of the linear problem in (27). In order to find all solutions of Eq. (26) in the domain Df 
the vector Fx k is repeatedly multiplied by the local root extractor 

ej(x,x^) 1 



1 — exp (— — x^H^) ' 

in which %W is the j-th solution of Eq. (26). In the repeated solutions of the transformed problem 



Fjx:= \Y\ej{x,x^>) )Fx = 0, J>1, (28) 




the process (27) is executed with a new Fx := Fjx. For every solution the process (27) is started 
many times with different x° and eo. Each time when J increases the derivatives /'(x^) are com- 
puted analytically and the matrices f' T ' (*Jt) ^ adaptively scaled [9]. The necessary last 
step of the method consists in a direct substitution check whether the Diophantine equation resid- 
ual vanishes exactly when the found solutions are rounded to integers. 

The system (24), (25) has been considered in Refs. [10] and [11] (pp. 285-286) as undecid- 
able. In Table 1 we present some examples in which the real-Diophantine system (24), (25) is 
solvable including the case when it is solvable over the primes (see cases 2, 3, 4 and 6). Here 
we shall describe in more detail two special examples which emphasize the crucial role of the 
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Table 1 : Real-Diophantine systems: when hi(xi) = Xi the solution is looked for among the integers 
while hi(xi) = p~ l (xi) means that the solution would be among the primes. 



case 


m 


n 


m 




source 


1 




3 


2 


Xi 


Pythagoras 


2 


xj+x^-xj-l = 


3 


2 


p- l (xd 


Sierpinski 


3 


x\ +X2+X3 +xf = n, neN 


4 


2 


p-'ixd 


Lagrange 


4 


9 , 

£ xf = n, n e N 

i=\ 


9 


2 


Xi, p~ l {xi) 


Waring, Khinchin 


5 


£ xf = n, neN 
i=\ 


19 


2 


Xi 


Waring 


6 


2 3 

(xi/X2) -{X3/X4) ~X 5 =0 


5 


2 


hi(xi) =p- l (xi), 
1=1,.. .,4, /z 5 (jc 5 )=x 5 


Fermat-Bache 



last step of the above method. In the first one we represent the prime number 5081 as a sum of 9 
cubes of primes (case 4 in Table 1 with hi{xi) = p~ l (xi)). We find two different solutions: 



2 x 2 3 + 3 x 3 3 + 5 3 + 2 x 11 3 + 13 3 , 
3 3 + 3 x 5 3 + 2 x 7 3 + 3 x 11- 



Notice that this problem has been solved as a real-Diophantine system on a machine with 16 
significant figures. The unknowns 2, 3, 5, 7, 11 and 13 in the first line of Eq. (29) have been 
found with 9 significant figures at residual | \f\ (jc) | |oo = 10~ 14 . A convergent process of the kind 
(27) has been build after 36 unsuccessful attempts which costed 5634 iterations with 4 different 
initial guesses ^0 combined with 9 different initial regularizators £q. The last step of the method 
yields an exact equality in Eq. (29). It would be interesting to investigate whether the number of 
primes which can be represented as the sum of 9 prime cubes is infinite. 
In the second example we consider the equation 

x\ \ I x 3 \ 



- =*s, (30) 

X2 J \X4 J 

where jc ; , i = 1,2,3,4 are sought as primes while JC5 as integer (special case of the unsolved 
Fermat-Bache problem in which the solutions are a rational pair (x\fx2, X3/X4) and an integer 
X5). The approximate solution found under conditions similar to those in the previous example, 
Eq. (29), 

xi =787.000011, x 2 = 348.99999357, jc 3 =457.00002128, 
x 4 = 1049.0000001, * 5 = 5.0024058062, 



leads to a nonzero residual after rounding to integers 

5 = 0.002405488240. 



787\ 2 / 457 \ 3 



349 J 



1049 
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This example is an illustration of the crucial importance of the last step of the method-the vector 
(787, 349, 457, 1049, 5) is not a true solution of the Fermat-Bache equation (30). 

The above method for the solution of Diophantine equations works because of a combina- 
tion of factors: autoregularization, SVD method, adaptive scaling, and because the solutions of 
Diophantine equations are well isolated. There is a semi-local convergence theory [6] in the non- 
degenerate case however no justification in the degenerate case is available by now. The methods 
of Refs. [4, 5] are applied to this problem with little success. 

Remark 2 The method (24), (25) for the solution of the Diophantine equation (24) does not 
contradict the negative solution of the 10-th Hilbert's problem [12] because our solutions are 
approximated in a bounded domain. 

4 The function p~ l (x) and the Riemann hypothesis 

Lemma 2 The functions p~ l (x) and k(x) are related by 

\p- l (x)-K(x)\<l VX>1. 

Proof: Let us assume that p n < x < p n+ i for some n. According to Eq. (20) the function p~ l (x) 
is strictly increasing, i.e., p~ l (x) < p~ 1 (p n+ \). On the other hand n(p n ) < n(x) < 7c(p n+ i) so 
that 

\p~ l (x) -k(x)\ < \p~ l {pn+\)-n{p n )\ = \n+l -n\ = 1, 

where we have used that p~ l (p n ) = n which follows from Theorem 1 and Eq. (5). 

□ 

Theorem 2 The asymptotics of p~ l (x) is the same as that for k(x) when x — > °°. 

Proof: Let us consider the relative difference of n(x) and Li(x). According to Lemma 2 we can 
write 

|K(*)-Li(*)| < \7t(x)-p-\x)\ | \p- l (x)-Li(x)\ < 1 | \p-\x)-U(x)\ 
Li(x) " Li(x) ' Li(x) — Li(x) Li(x) 

Therefore in the limit x^«we can ignore the term 1/Li(jc) and investigate p~ l (x) instead of 

k(x). 

□ 

4.1 Differential equation and the von Koch estimate 

Let us consider the following function 

£ ^M_u( f ) x> 

wxm(x) 
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where Li(x) is defined in Eq. (1). According to the von Koch estimate (see [1], pp. 90) the 
Riemann hypothesis is equivalent to the statement that K(x) is asymptotically constant, i.e., 

lim K{x) = const. -<=>- RH. 



Because the function p 1 (x) is continuously differentiable we can write the following differential 
equation for K 

v ; \2x xln(x) / v ; v^ln(x) 

The derivative dp~ l (x)/dx is strongly oscillating as shown on Fig. 2, however it is restricted 
between and 1 according to Eq. (20). Therefore we shall consider the solution of Eq. (32) in the 
interval p n <x< (p n + p n+ \)/2, n ^ <*>, and shall use the fact that (see Eqs. (14) and (16)) 



dp (x) dp (x) 



dx 



dx 

■* — Pn 



= ^J—t- (33) 

X=(pn+Pn+l)/2 l Pn+l 



Thus, for x — > p„ we can substitute dp /dx = 1 in Eq. (32), neglect the term 1 / ln(x) in the limit 
x — > °° and solve the equation 

x 'W = -(^ + JTSw) x W + ^W' (34) 

The general solution of this equation can be written as 

K(x) = -y^+- c ' n 



ln(x) ^xln^)' 

where the first term in the right-hand- side is a partial solution of the inhomogeneous equation (34) 
while the second one is the general solution of the homogeneous equation and c' n is a constant. 

At the right-hand border* — > (p n + p n+ \)/2, we can neglect, for x — > °°, the term dp~ l /dx = 
(2p n + 1) _1 and keep only l/ln(x) assuming that p„ ~ with a > 0. In this case we should 
solve the equation 

Tx + x^) K{x)+ ^/xh?( x y 

The general solution of (35) is again the sum of a partial solution (the first term bellow) of the 
inhomogeneous equation and the general solution (the second term) of the homogeneous one 

/ s Li fx) c" n 

K(x) = =ff- + ■ 



K ' to = -(- + J *to + "FriTT- (35) 



v^cln(x) v^cln(x)' 

Substituting the logarithmic integral with its leading term for x — > °°, i.e., Lifx) ~ %/ lnfx) we can 
finally write 



K{x) 



V|_I_C%_ Ph+Ph+1 



It is tempting to regard the general solution of the homogeneous equation as subleading in the 
limit x — > oo and the first terms as expressing the oscillations of K(x) close to the borders of the 
considered intervals. However, let us note that the constants c' n and c" n might depend on the 
primes gap p n which on its own depends on p n and this last dependence is currently unknown. 
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4.2 The l'Hospital rule 

Here we shall consider the limit 

Iim p-Hx)-Li(x) = 1 . m p^(x)_ 1 
Li(jc) x-t- 00 Li(jc) 

Because the function p^ 1 (x) is differentiable we can apply the l'Hospital rule if the limit 

lim = i im i n ( x ) ap W 

ln(x) 

exists. Now let us show that if the RH is true then this limit does not exist. Indeed, let us choose 
two subsequences of x — > °°, namely 

(/) x = p m n ->• oo, 

(u) x= , n^oo. (36) 

Then, using again the values (33) of the derivative over (i) and (ii) we get 

(0 ln(*)^^~ln(p„)^oo, 

lnW ^ ^TT^°- (37) 

i +e 

The second limit follows from the statement that if the RH is true then p n ~ p% for any £ > 
when n — > oo [2]. Most of the current estimates of the primes gap p n lead to the non-applicability 
of the l'Hospital rule. Nevertheless we cannot be sure until a rigorous estimate is found. 

5 Conclusions 

We have constructed a pair of diffeomorphisms p(x) and p~ l (x) which interpolate the prime 
series and the prime counting function, respectively, which are convenient for both numerical and 
analytical applications. To the best of our knowledge this is the first differentiable and invertible 
interpolation of the prime series. 

The function p~ l (x) can be effectively used for the solution of Diophantine equations which 
can be exploited in many cases where the other methods do not work and could be particularly 
useful when Diophantine equations are subsystems of more complex real systems. 

Because p~ l (x) has the same behavior as k(x) for x — > °o it could give more information about 
the asymptotic and non- asymptotic distribution of primes. Perhaps, this could be used to draw 
some conclusions about the Riemann hypothesis when more information about the primes gaps 
becomes available. 
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